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We report tilted-field magnetotransport measurements 
of two-dimensional electron systems in a 200 A-wide 
Alo.13Gao.87As quantum well. We extract the energy gap for 
the quantum Hall state at Landau level filling = 1 as a 
function of the tilt angle. The relatively small effective Lande 
g-factor {g ~ 0.043) of the structure leads to skyrmionic exci- 
tations composed of the largest number of spins yet reported 
(s ~ 50). Although consistent with the skyrmion size ob- 
served, Hartree-Fock calculations, even after corrections, sig- 
nificantly overestimate the energy gaps over the entire range 
of our data. 
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In two-dimensional electron systems (2DESs), the 
quantum Hall effect (QHE) at Landau level filJfijtg fac- 
tor v ~ has attracted much theoreticalETQ and 
experimentaflilj attention. At this filling, the Coulomb 
(exchange) energy is so influential that the QHE excita- 
tion gap is more than an order of magnitude larger than 
the single-particle Zeeman energy, the gap expected if 
the Coulomb energy were "turned off." In fact, given a 
small enough Zeeman energy, the interplay between these 
two energies leads to a lowest-lying charged excitation, 
called a skyrmion, composed of electrons arranged in a 
canted, nearly-parallel spin-texture. Properties of this 
excitation, such as its energy gap and physical extent, 
are determined by the ratio g = \g\fJ,BBtot/{e'^ /^f-s) of 
the single-particle Zeeman energy, which limits the num- 
ber of spin-flips in an excitation, to the Coulomb energy 
which favors local ferromagnetic ordering, (e is the di- 
electric constant, £3 — \/tic/eBj_ is the magnetic length, 
g is the effective Lande g-factor, iib the Bohr magneton, 
and Btot and B± are the total applied magnetic field and 
the component perpendicular to the layer plane respec- 
tively.) 

The limit of g ^ is of particular interest where, 
in an ideal system, the excitation gap is predicted to 
exist even in the absence of Zeeman energy and the 
skyrmion size diverges — s — » 00. (We choose s to de- 
note the total spin of a thermally-activated skyrmion- 
antiskyrmion pair; the effective spin of a single skyrmion 



or antiskyrmion would therefore be s/2 if particle sym- 
metry holds.) Experiments employing diverse techniques, 
such as optically-raiHiped nuclear magnetic resonance,Ll 
magnetotransport ,eHU and magneto-optical absorption 
spectroscopya in GaAs 2DESs, where \g\ ~ 0.44, have 
yielded s 7 to 9. By using hydrostatic pressure to 
tune g, Refs. p]-[T^ can access the regime g ^ where 
they extract a larger number of spin-fiips (s = 36) from 
the temperature dependence of their magnetotransport 
data. Unfortunately, this technique requires a separate 
cooldown for each g which leads to a different disorder 
potential every time. Furthermore, since applying hy- 
drostatic pressure lowers the density of the 2DES, to 
compensate, the sample must be illuminated. Because 
of these complications, a controllable, systematic study 
using pressure is non-trivial. 

In this paper, we report the observation of the largest 
skyrmions yet reported (s ~ 50) by using an alternate ap- 
proach to access the low g regime. In bulk Al^rGai-^^As, g 
increases monotonically from g = —0.44 at a; = (GaAs) 
to g w +0.5 at X = 0.35 vanishing at a; ~ 0.13.eI We 
fabricated a wafer with a 200 A Alo.13Gao.87 As quantum 
well bounded first by a thin 12.6 A AlAs layer followed 
by thick Alo.35Gao.65As barriers on each side. Grown 
by molecular beam epitaxy, this symmetric structure is 
modulation doped with -Si. While many experimental 
techniques measure (;,t3 l2I determining an extremely low 
g is difficult and subject to great relative—Lmcertainty. 
From calculations utilizing the Kane modelj£j'E3 our best 
estimate of the g-factor for this wafer is g — 0.043±0.01CL 
an order of magnitude lower than that in bulk GaAs.o 
The slightly positive value of g is attributed to the spill- 
over of the electron wavefunction into the barrier re- 
gion where g is positpjSy-as well as the non-parabolicity 
of the energy bands.cJ'EJ In the right inset to Fig. 1, 
we show a calculation of the longitudinal g-factor for a 
symmetrically-distributed electron system as a function 
of Al concentration a; in a 200 A Al^jGai-^jAs quantum 
well with barriers as specified above. We have taken the 
g-factor to be isotropic; we will address this assumption 
later. 

We measure two samples (identified as A and B) 
from different parts of the wafer with mobility /i « 
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5 X 10"* cm^/Vs in a Van der Pauw geometry. Sam- 
ples A and B have total areal densities n = 1.37 x 10^^ 
and 1.28 x 10^^ cm~^, respectively. We collect the low- 
temperature magnetotransport data in a dilution fridge 
and a ^He system. In our experiment, we extract the 
excitation energy gap (Ai) for the QHE at filling v = 1 
from the temperature-dependence of its longitudinal re- 
sistance {Rxx) minimum. We gather Ai for several g 
by tilting the sample and thus applying a magnetic field 
-Btot at an angle 6 with respect to the normal of the 
sample plane. This technique allows the Zeeman energy 
(oc i?tot) to be tuned in situ while tke other parameters 
in the system are nearly unaffected.u Unlike the experi- 
ments in Refs. |l^-|lj, the areal density n, the disorder, 
and the Coulomb energy e? jeiB remain constant for dif- 
ferent values of g leading to a relatively straightforward 
analysis of our data. By finding the energy gap for sev- 
eral angles Q (or equivalently 5), we can determine the 
number of spin-flips involved in an excitation since any 
change in the gap is almost entirely attributable to the 
change in the Zeeman energy contribution. As in Ref. ^ 
and Refs. pl|-p^, we use the formula s — dAi/dg, where 
Ai = Ai/(e^/e£s) is Ai normalized by the Coulomb en- 
ergy, to extract the number of spin-flips in an excitation. 

The traces in Fig. 1 attest to the high quality of the 
sample as well as the small value of the g-factor. The 
longitudinal magnetoresistance of sample A is plotted 
for temperatures T = 20 mK and 735 niK. The 20 mK 
trace exhibits minima for QHE states with even inte- 
ger fillings as high as = 42 at B± = 0.13 T mark- 
ing the lower bound for the Shubnikov-de Haas (SdH) 
oscillation regime. (See the left inset to Fig. 1.) We 
believe that impurity scattering, rather than alloy scat- 
tering, is the dominant mechanism limiting the mobil- 
ity in our AlGa A s. quantum wells. A Born approxima- 
tion treatmentOcS estimates the alloy scattering mobil- 
ity limit to be « 4 X 10^ cm^/Vs, an order of magnitude 
larger than our measured mobility. Impurity scattering, 
however, may explain the low measured mobility. It is 
known that AI, a relatively reactive element, incorpo- 
rates impurities in AlGaAs layers during growth.LlI Such 
a mechanism is consistent with the relatively low mo- 
bility (even after considering the occupatipM of multiple 
ellipsoids) observed in AlAs quantum wcUsES where there 
is no alloy scattering. 

Although minima corresponding to QHE states exist 
for many even integer fillings, there are none for odd in- 
tegers other than v — 1 and 3. At higher odd fillings, 
the influence of the Coulomb exchange energy is progres- 
sively diminished since the states occur atJower and 
the fraction of electrons affected is 1/i^;l3 instead, the 
excitation gap at higher odd v is determined primarily 
by the competition of the single-particle Zeeman energy 
and disorder-broadening of Landau levels. In our system 
the Coulomb exchange energy appears to be significant 
enough to overcome disorder-broadening only for the odd 
integer fillings v = \ and 3. In fact, calculations which 



consider the finite width of the 2DESl3 predict skyrmionic 
excitations for our sample at both these fillings. Unfortu- 
nately, the measured excitation gap at = 3 is only 
K which is of the order of the Landau level broadeningEEl 
making further analysis inconclusive. 

We now discuss our data ioY v — 1 . Shown in Fig. 2 are 
the Arrhenius plots of Rxx minima for the v = \ QHE 
in sample A for three angles in the temperature range 1 
< T < 4 K. We extract the activation energy Ai from 
the slope of a best-fit line (dashed line) to the data using 
the relation Rxx ~ exp(-Ai/2r). For 6* = 0, 49.7, and 
67.6°, we have Ai = 15, 20, and 25 K respectively. Our 
measured Ai is a monotonically increasing function of 9 
lying in the range 13 < Ai < 25 K as shown in the inset 
to Fig. 2. Since the data for both samples A and B are 
qualitatively very similar, we focus on sample A below. 

In Fig. 3, we plot Ai = Ai/(e^/e^B) vs g and discuss 
our results in light of other experiments and theoretical 
calculations. The data from sample A, which occupies 
the extreme lower left portion of the figure, is expanded 
in the inset to Fig. 3. An asymptote (dashed line) fit 
to the lower range of the data reveals s = 50.2 ± 1.0. 
Similar analysis on sample B yields s = 49.2 ± 2.1. Com- 
pared to our experiment, Ref. |^ explores higher g in GaAs 
samples represented by various closed symbols in Fig. 
3. In the case of the pressure-tuned data (not shown) 
from Refs. [Tl|-|T^, g is in the same vicinity as our data, 
although there is more scatter in the reported Ai, pre- 
sumably because of variations in disorder as discussed 
before. A noteworthy commonality in these experiments 
is the range of the measured skyrmionic excitation gap 
(Ai < 0.33). 

We now contrast our experimentally obtained Ai with 
calculations. The top solid curve in Eiff. 3 represents the 
results of Hartree-Fock calculationaa'El for the skyrmion 
excitation energy gap in an ideal, infinitely-thin 2DES, 
i.e. one with an electron probability density width w ^ Q. 
The calculated skyrmion gap declines steeply as 5 — » 0, 
reflecting the decreasing cost in Coulomb energy for an 
excitation with an increasing degree of nearly-parallel 
spins. On the other hand, the (exchange-enhanced) sin- 
gle spin-flip excitation gap (dashed line) expected in the 
absence of skyrmions has a constant slope corresponding 
to s = 1. At g = 0, the skyrmion gap is ^-y/f^e^/e^s, half 
the single spin-flip excitation gap. The skyrmion remains 
the favored excitation for g < gc = 0.054 (marked by a 
vertical arrow). 

Note the striking discrepancy between the calculations 
and the experiments. The calculated skyrmion gap for 
the ideal case is a factor of 4.3 to 6.7 larger than the 
experimental data! The ideal case, however, ignores im- 
portant effects such as finite thickness correction (FTC), 
Landau level mixing (LLM), and disorder-broadening of 
Landau levels — all of which reduce the energy gap. Un- 
fortunately, since no calculation currently treats these 
three corrections simutaneously for skyrmions, we must 
consider them in cumulative succession. Figure 3 in- 
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eludes Hartree-Fock calculations with FTCH for layer 
thickness w = 0.43 £5, appropriate for our sample. 
(We determine w by fitting a gaussian function to the 
electron probability density from a self-consistent local 
density approximation calculation.) By softening short- 
range interactions, the FTC reduces the predicted gap 
by w 30%. To assess the effect of LLM, we first focus on 
the exchange-enhanced single spin-flip excitation gap cor- 
rected for FTC and LLM. Existing calculations are not 
in quantitative agreement; based on the trends in Ref. 
pT| and Ref. we estimate the corrected single spin-flip 
gap to be about 0.58 and 0.69 e'^/eis respectively for our 
sample (plus the Zeeman energy, Iffl/isBtot)- And, if the 
role of disorder is limited to the disorder-broadening of 
Landau-levels, we expeqt-|the predicted gap to diminish 
by only F ~ 0.06 e^/eisS^ We then deduce the gap in the 
large skyrmion limit (small g) by shifting the w — 0.43 £3 
curve for the skyrmion excitation gap by a constant to 
match (at g > cjc) the corrected single spin-flip gaps es- 
timated. The shifted curve still overestimates our exper- 
imental gaps for = 1 by a factor of 1.3 to 1.5 and 1.8 
to 2.3 for Ref. |3^ and Ref. ^ respectively. 

Although the absolute values of the calculation for the 

= 1 QHE cannot be reconciled with the experimental 
data, the size of skyrmions (from the slope of the curve) 
predicted by the calculation agrees with our data. This 
agreement is evidenced by the w — 0.43 £3 curve shifted 
down by 0.46 e^jels (shown by the dotted line in the 
inset to Fig. 3) which fits the entire lower range of our 
data remarkably well. In fact, we can use the calculation 
to check the validity of our value for the g-f actor. We find 
that this agreement is valid only in a very narrow range 
of assumed_.values for g which includes our estimate of 
g ~ 0.043.e3 This congruity may be viewed, perhaps, as 
an independent confirmation of g-factor in our sample. 

Thus far, we have interpreted our data assuming an 
isotropic g-factor for our sample. In general, however, the 
(^-factor can be anisotropic in confined systems, with gg 
and gt denoting the longitudinal and transverse compo- 
nents of the g-factor with respect to the growth axis (for 
a review, see Ref. |2| and references therein) . The elec- 
tron g-factor anisotropy is governed by the low-symmetry 
electron quantum confinement and changes strongly with 
the quantum well width; it can be qualitatively estimated 
from the energy splitting between the light- and heavy- 
hole bands. Based on a time-resolved photoluminesceneej 
experiment measuring electron spin quantum beats,Lil 
Le Jeune ei al. conclude that the g-factor is indeed 
anisotropic for narrow GaAs quantum wells bounded by 
Alo.3oGao.7oAs barriers, as found in Refs. |20| , |35| ; however, 
in quantum wells 120 A and wider, the anisotropy van- 
ishes. We note here that the g-factor anisotropy may be 
reduced as the electron kinetic energy in the 2DES, the 
thermal energy, or localization energy (because of imper- 
fections or magnetic field) become comparable in value 
to the quantum confinement energy. 

A simple Kane-model calculation for electrons at 
the bottom of the first subband in our system yields 



a transverse component of the 51-factor gt = 0.085. 
And, if we reinterpret our data so that g[6) — 

\J gjcos^ {6) -I- 5^sin^ {9) , then s from the asymptote to 
the lower range of the data is reduced — s = 19. How- 
ever, experimental findings mentioned above support an 
isotropic (^-factor for the parameters in our sample. Our 
system with a quantum well of width 200 A and a finite 
2DES concentration should exhibit even less of a ten- 
dency toward (7-factor anisotropy. And, as already noted, 
the size of skyrmions from the calculations and our data 
are no longer consistent if the anisotropic g-factor is as- 
sumed. Therefore, we believe that the isotropic value of 
0.043 is the best estimate for the g-factor in our samples. 

In summary, we have focused on the thermal excitation 
energy for v — 1 gathered for several tilt-angles from 
magnetotransport measurements of two-dimensional 
electron systems in a 200 A-wide Alo.13Gao.87As quan- 
tum well. In this structure with a small g-factor [g ~ 
0.043), we observe skyrmions of the largest size yet re- 
ported (s ~ 50). The magnitude of the energy gaps 
measured are consistent with those from other experi- 
ments. And, while matching the experimentally deter- 
mined size of skyrmions, Hartree-Fock calculations, even 
after treatment for corrections, significantly overestimate 
the energy gaps in our data. Understanding this dispar- 
ity requires further studies in the large skyrmion regime. 
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FIG. 3. Normalized activation energy Ai — Ai / {e^ /sIb) 
vs normalized Zeeman energy g = (;/is-Btot/(e^/e£s) from ex- 
periments and calculations. The experimental data are from 
sample A (+) and Ref. 8 (closed symbols). In the inset, the 
asymptote (dashed line) fit to the lower range of sample A 
data reveals s = 50.2 ± 1.0. The results of Hartree-Fock cal- 
culations for a 2DES with zero layer-thickness {w — 0) and for 
w = 0.43 is are also shown in the main figure. In the inset, 
the w = 0.43 £b skyrmion excitation gap (dotted line) shifted 
down by 0.46e^/e£B matches the lower range of sample A 
data. (See text for details.) 



FIG. 1. Magnetoresistance traces for sample A at T = 20 
and 735 mK which show QHE states for only i/ = 1 and 
3 among odd f. The left inset shows the onset of Shub- 
nikov-de Haas oscillation at v = 42 {B± = 0.13 T). The 
right inset shows a calculation of the longitudinal g-factor 
c/£ for a 200 A AUGai-^jAs well bounded by Alo.35Gao.65As; 
gi = 0.043 ± 0.010 for x = 0.13. 



FIG. 2. Arrlienius plot of Rxx minima for the i> = 1 QHE 
state at several angles 6 for sample A. The activation energy 
Al is plotted in the inset vs 9 for both samples A and B. 
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Fig. 1 "Large Skyrmions in an Al^ j^Gagg^As quantum well" 
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Fig. 2 "Large Skymions in an Al Ga As quantum well" 
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Fig. 3 "Large Skyrmions in an Al Ga As quantum well" 




